Neutrino mixing as a source for cosmological constant 
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We report on recent results showing that neutrino mixing may lead to a non-zero contribution 
to the cosmological constant. This contribution is of a completely different nature with respect to 
the usual one by a massive spinor field. We also study the problem of field mixing in Quantum 
Field Theory in curved space-time, for the case of a scalar field in the Friedmann-Robertson- Walker 
metric. 
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I. INTRODUCTION 

The problem of cosmological constant is currently 
one of the most challenging open issues in theoretical 
physics and cosmology. The main difficulty comes from 
the fact that, when estimating its value as a gravita- 
tional effect of vacuum energy, the numerical results are 
in strong disagreement with the accepted upper bound 
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cm 



A < 10 

We show here that a new contribution to the vacuum 
energy and therefore to the cosmological constant may 
arise from neutrino mixing . The contribution we find 
comes from the specific nature of the field mixing and 
is therefore of completely different origin with respect to 
the ordinary vacuum energy contribution of a massive 
spinor field. 

Indeed, it has been shown 0-01 ^^^^ the vacuum for 
fields with definite masses is not invariant under the field 
mixing transformations and in the infinite volume limit, 
it is unitarily inequivalent to the vacuum for the fields 
with definite flavor. In the case of neutrinos, this re- 
sults in a condensate of neutrino-antineutrino pairs, with 
a density related to the mixing angle (s) and mass differ- 
ence(s) among the different generations. Phenomenolog- 
ical consequences of such a non-trivial condensate struc- 
ture of the flavor vacuum have been studied for neutrino 
oscillations |1] and for beta decay 

In this paper we consider the case of two flavors and 
Dirac neutrino fields, although the conclusions we reach 
can be easily extended to the case of three flavors and 
Majorana neutrinos [Tllll^ . 

We also include a preliminary study of mixed (bosonic) 
fields in a curved background, for the case of FRW metric. 

In Section II, we shortly summarize the main results 
for neutrino mixing in QFT. In Section III, we compute 
the neutrino contribution to the cosmological constant 



and estimate its value by using the natural scale of neu- 
trino mixing as a cut-off. The result turns out to be 
compatible with the currently accepted upper bound on 
A. In Section IV, quantum fields and mixing relations are 
analyzed in expanding universe. Section V is devoted to 
conclusions. 



II. NEUTRINO MIXING IN QUANTUM FIELD 
THEORY 

The main features of the QFT formalism for the neu- 
trino mixing are summarized below. For a detailed re- 
view see [l3|. For sake of simplicity, we consider the 
two flavor case and we use Dirac neutrino fields. The 
Lagrangian density describing the Dirac neutrino fields 
with a mixed mass term is: 



(1) 



The relation between Dirac fields ^'/(x), eigenstates of 
flavor, and Dirac fields '^m{x), eigenstates of mass, is 
given by 



U is the mixing matrix 



(2) 



l(=\ ''''' (3) 
sm cos (> ' ^ ' 



being 9 the mixing angle. Using Eq.(|3Jl, we diagonalize 
the quadratic form Eq.lQJ, which then reduces to the 
Lagrangian for the Dirac fields ^'^(x), with masses m^, 
i = 1,2 : 



C{x) = *„(a;) {i 0- Md) *™(a;) , 



(4) 



where — diag{mi^m2). The mixing transformation 
^ can be written as 3] 
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y,{x) = G,\t)v,{x) Ge{t), 



(5) 
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where (cr, i) 
given by 



(e,l),(^,2), and the generator Gg{t) is 



Geit) =exp e J d^x(^i'lix)Mx)-4{x)Mx)) ■ (6) 

The free fields i>i (i=l,2) are given, in the usual way, 
in terms of creation and annihilation operators (we use 
t = xo): 



with 1 = 1,2, 



(2^)1 



-k,i 



«k..(i) 



e I'k,; 



(7) 
(8) 



and ujk,i = ^/k^ + mf . 

The mass eigenstate vacuum is denoted by |0)m: 
"^^110)771 = /3ki|0)m = 0. The anticommutation rela- 
tions, the wave function orthonormality and complete- 
ness relations are the usual ones (cf. Ref. 0]). 

The flavor fields are obtained from Eq. lO: 



E 



k,o- 



"k,»(i)"k.<T(i) 



(9) 



with (cr, i) = (e, 1), (/i, 2), with the flavor annihilation 
operators deflned as 

Ge\t) «M Ge{t) 



^k,<T 



it) 



(10) 

/3!:k,.(i) ^ G,\t) (3%Ge{t). (11) 
They annihilate the flavor vacuum |0(t))/ given by 

W))f ^ G^\t) |0)„, . (12) 

In the infinite volume limit, the vacuum |0(^)) / for 
the flavor flelds and the vacuum \0)m for the flelds with 
definite masses are unitarily inequivalent vacua [H, 3- 

One further remark is that the use of the vacuum state 
|0),jj in the computation of the two point Green's func- 
tions leads to the violation of the probability conservation 
0. The correct result is instead obtained by the use of 
the fiavor vacuum \0)f, which is therefore the relevant 
vacuum to be used in the computation of the oscillation 
effects. We will thus use |0)^ in our computations in the 
following. 

The explicit expressions for the flavor annihila- 
tion/creation operators in the reference frame k = 
(0,0, |k|) are Q: 



^k,2 



e'-Fk(t)n,2 



Q^k,e(^) — cos0al^ I + sin (u^(t)a'^ 
"k.^W = cos 61a;, 2- 

P-Keit) = cose f3L^^, 

P-K^t) = cos 0/3:^,2 - sine (c/k(0/3!:k,i + ^''y^{t)<^Zi) 

(13) 



sme (U^{t)al^^~ €'V^{t)l3\^ 
+ sin^^({/^(^)/3:k,2-e'-t/k(^)a'■t 



k.2 




Lo£(|k| 



FIG. 1: Fermion condensation density |Vk|^ as a function of 
|k| and for sample values of the parameters mi and 7712. Solid 
line: mi = 1, m2 = 100; Long dashed line: mi = 10, m2 = 
100; Short dashed line: mi = 10, m2 = 1000. 



where C/k and Vk are Bogoliubov coefficients given by: 
14(t) = l^kl e*('^'''=+"'=-i)*, (14) 



|f/k 



/ LUkS + mi \ 2 / U}k^2 -I- TO2 V 



(15) 



X 1 



IkP 



m = 



+mi\^ ( uok,2 + m2 



Ikl 



2w*:,2 
Ikl 



with 



{U!k,2+m2) {uJk,! + mi) J ' 



|C/kP + |VkP = 1. 



(16) 



(17) 



The function |Vk| is related to the condensate content 
of the flavor vacuum |3| as: 

e,M(0|<.aM|0)e,M - sin^e 1141' , 1,2 (18) 

with the same result for antiparticles. 

In Fig.l, we plot the fermion condensation density 
|Vkp in function of |k| and for sample values of the pa- 
rameters mi and 1712- |Vkp is zero for mi = m2, it has 
a maximum at |k| — ^m^m^ and, for |k| 3> ^Jm^ra^^ it 
goes like \Vv\^ ~ (m, - mJV(4|k|2). 

From the study of the current algebra, we obtain the 
fiavor charge operators that, in terms of flavor operators, 
are: 

QAt) = E {<At)<At) ~ ii%,At)PWM) , (19) 

k,r 
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with (7 = e, /i. 

At time t — 0, the vacuum state is |0)e,/j and the one 
electron neutrino state is defined as: 

We) = aL!e|0)e,A' 

Considering the flavor charge operators, defined as in 
Eq. p9|l . We then have (in the Heisenberg representation) 



where G is the gravitational constant. 

The energy-momentum tensor density T^j^ is obtained 
by varying the action with respect to the metric g^i/i 



T — 



8S 



.(o|ge(i)|o) 



e,M-e,M (0|QMW|0>e,M=0, (20) 



where the action is 



~gC{x)d'^x. 



(27) 



(28) 



(21) 



{«Ic.eW,<e(0)}f +|{/3:keW><e(0)} 



(22) 



= {«L,,(i),<e(0)}f +|{n,,W,<e(0)} 

Charge conservation is obviously ensured at any time: 

QleW + Qljt) = 1- (23) 

The oscillation formula for the flavor charges are then 



Ql^it) = l-sin^(20) \U^f sin' Z'^'-' t 



,2.2/ ^k,2 +(^k,l, 

|Vk| sm I 1 



QtJt) = sin^(20) \U^f sin' ( ^^M_^t 



Vk sm I z 1 



(24) 



(25) 



This result is exact. There are two differences with re- 
spect to the usual formula for neutrino oscillations: the 
amplitudes are energy dependent, and there is an addi- 
tional oscillating term. 

2 2 

For |k| 3> ^/mim2 we have \Uk\ — > 1 and |Vk| — > 

and the traditional formula is recovered. 

Similar results are obtained in the case of boson fields 
[To|. For the rj — r]' system, the correction may be as large 
as 20%. 



III. NEUTRINO MIXING CONTRIBUTION TO 
THE COSMOLOGICAL CONSTANT 

The connection between the vacuum energy density 
(Pvac) and the cosmological constant A is provided by 
the well known relation 



{Pvac 



A 

4^' 



(26) 



In the present case, the energy momentum tensor den- 
sity is given by 

(^m(x)7^(x)'i^^^'™(a;)) , (29) 
where I? is the covariant derivative: 



(30) 

being 7"^ the standard Dirac matrices, and ^' D i.'i' — 
^D'i> — (D^)'^. Let us consider the Minkowski metric, 
we have 



oo(a;) = ^ : (^"S ^{x)jo^ o"^ mi^)) : (31) 



where : . . . : denotes the customary normal ordering with 
respect to the mass vacuum in the flat space-time. 

In terms of the annihilation and creation operators of 
fields i>i and 1^2 : the energy- momentum tensor 



Too - J d^xTooix) (32) 

is given by 

= E / ^'k^fc. + n/^k,) , (33) 

with i = 1,2. 

Note that Tqq is time independent. 

The expectation value of Tqq in the flavor vacuum |0) / 
gives the contribution (p^ac) of the neutrino mixing to 
the vacuum energy density is: 



/(0|E^oo(0)|0)/ = (/^"- 



(34) 



Within the QFT formalism for neutrino mixing, we 
have 

f{0\T^l^\0)f = fm\T^l^m)f (35) 
for any t. We then obtain 

i i.r 

+ fmZPlMf)- (36) 



4 



By using Ea. (fT5|) . we get 

/(0|I]T^oo^(0)|0)/ = SsinH J d^k{uk,i+iOk,2)\Vk\^ 

i 

- {pT^Dvoo, (37) 

i.e. 

{pZc) = 327r2sin2 / dkk\ujk,i + ^^^,2)^^ (38) 
Jo 

where the cut-off K has been introduced. Ea. (|38|l is our 
result: it shows that the cosmological constant gets a 
non-zero contribution induced from the neutrino mixing 
0. Notice that such a contribution is indeed zero in 
the no-mixing Hmit when the mixing angle 9 — Q and/or 
TOi = m2. Moreover, the contribution is absent in the 
traditional phenomenological (Pontecorvo) mixing treat- 
ment. 

We may try to estimate the neutrino mixing contribu- 
tion by making our choice for the cut-off. If we choose 
the cut-off proportional to the natural scale appearing in 
the mixing phenomenon ko ~ ■^mim2 3]: using K ^ fcoi 
mi = 7 X lO-^eT^, ma = 5 x IQ-'^eV, and sin^ 9 ~ 0.3 
[3 in Eq.lEHI), we obtain 

= 0.43 X IQ-^^GeV^ (39) 

Using Ea. H26|l . we are in agreement with the upper 
bound for A: 

A - IQ-'^^cm-^, (40) 

Another possible choice is to use the electro-weak scale 
cut-off: K w lOOGeF. We then have 

ifC^) = 1.5 X lO-'^GeV' (41) 

and 

A - 10-24cm-^ (42) 

which is, however, beyond the accepted upper bound. 

In a recent paper [13 , it was suggested the cut-off scale 
given by the sum of the two neutrino masses, K — mi -|- 
m2. 

For hierarchical neutrino models, for which m2 S> mi, 
we have, in this case, K 3> y^mimj, and thus, if we as- 
sume that the modes near the cut-off contribute mainly 
to the vacuum energy, and take into account the asymp- 
totic properties of Vk: 

l^^|,^(m2_W k»V^2 (43) 

we obtain: 

{P7ac) STTsin^ 6'(m2 - mi)2(m, + m x 



and then 

(p-i-)«sin2 0(Am2)2 (45) 

in the limit m2 2> mi. 

In Ref . TE\ , the corresponding Am^ is given by the so- 
lar neutrino data: Am^ ~ 10~^eV^, resulting in a con- 
tribution of the right order. 

IV. QUANTUM FIELDS AND MIXING IN 
EXPANDING UNIVERSE 

In this Section we present a preliminary study of mixed 
fields on a time-dependent gravitational background. 
For simplicity, we consider the case of neutral scalar 
fields, the case of fermionic fields and neutrino oscilla- 
tions will be considered elsewhere. 



A. Free fields in expanding universe 

We start by quantizing a free neutral scalar field 
in the Friedmann-Robertson- Walker (FRW) space-time 
with flat spatial sections, characterized by metrics of the 
form 

ds^ = gi^^dx^'dx" = dt^ - a2(^)^^x^ (46) 

where a{t) is the scale factor. 

A flat FRW space-time is a conformally flat space-time. 
Indeed, by replacing the coordinate t by the conformal 
time r/, 




where to is an arbitrary constant, the line clement takes 
the form 

ds^ = a^{'n)[dri^ ~dyi% (48) 

where a{ri) is the scale factor expressed through the new 
variable rj. Introducing the auxiliary field x = o,{r])(j), it 
is possible to show that the evolution of a scalar field 
in a flat FRW metric is mathematically equivalent to 
the dynamics of the auxiliary field x i"^ the Minkowski 
metric ,1(|. The information about the influence of grav- 
itational field on (j) is contained in the time-dependent 
mass meff{r]) defined by 

a" 

tnlffivi) ^ m^a^ : (49) 

" a 

where the prime ' denotes the derivative with respect to 
V- 

The field x can be quantized in the standard fashion 
by introducing the equal time commutation relations 

[X(x,77),7r(y,77)] =*(53(x-y), (50) 
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where tt = x' is the canonical momentum. 
The Hamiltonian of the quantum field x is 



H{r,) = ^ I cP^ [tt^ + (Vx)' + m%f{T^)x^] . (51) 



Note that the energy of the field x is not conserved; this 
leads to the possibility of particle creation in the vacuum. 
The energy for new particles is supplied by gravitational 
field. 

The field x is expanded as 



X(x,?7) 



(vl{ri)a^ + Vk{'q)al^ , (52) 



(27r)i a/2 

where the mode functions Vkirf) obey the equations 



'"k + ^k{v)vk = 0, ujkiv) = \/k'^ + mlff{r)) (53) 
and satisfy the following normalization condition 

v'kVl - Vkv'C = 2z. (54) 

We have Vk{ri) — v^k{^) as follows from the relation 
(Xk)* = X-k and (ck)* = aj^, where 



Xk(ry) = ^ (t'fe('?)ak + Vk{il)a}_C) 



(55) 



The annihilation operators are given in terms of Xkl??) 
and Vk{vi) by: 



flk 



= V~2 



<Xk - t'fcxL 



Note that at is time-independent; the commutation re- 
lations are: 

[«k, al] = 53(k _ k'), K, Ok'] = [4, 4'] = 0- (57) 

The vacuum state annihilated by ak is denoted by 
|0(ry)): ak|0(77)) = 0. 



B. The instantaneous vacuum 

The Hamiltonian Ea. (|51|l is time dependent, then we 
can define an "instantaneous" vacuum by selecting an ar- 
bitrary time 7^0 and defining the vacuum |0(77o)) as the 
lowest energy eigenstate of the Hamiltonian H{r]o) com- 
puted at the time rjo- 

The mode functions w/c(?7o)i corresponding to the vac- 
uum 10(770)), are obtained by computing the expectation 
value {0{ri)\H {r]o)\0{r])) in the vacuum state |0(ry)) de- 
termined by arbitrary chosen mode functions Vk(j]) and 
then minimizing that expectation value with respect to 
all possible choice of Vk{r]) fl^ . 



The Hamiltonian Ea. (|51(l expressed through the an- 
nihilations operators at defined by the arbitrary mode 
functions Vkir]) is 



v'k'iri)+u^liv)vl{ri))aial^ 



+ 



( 



(58) 



Since ak|0(ry)) = 0, we have 

+ u;l{rj)\v,irj)\A 

and the density energy is 

p^\[ d'k(Kirjr+u;lirj)\v,irjr) . (59) 

At fixed value of the momentum fc, if 1^^(779) > 0, it 
is possible to show that the mode functions Vk{r]) that 
minimize p satisfy the following initial conditions at 77 = 
770 .16]: 



Vk{r]o) 



1 



^kim) = iy/uJkM- (60) 



V^kivo) 

The mode functions satisfying the Ea. H6n|l define 
the annihilation operators of the vacuum |0(77o)): 
0^10(770)) = through which the Hamiltonian H{r]o) at 
time 770 is expressed as 



(56) Hifjo)= J d'kiJkim){4^< + ^S'{0)). (61) 



The zero point energy density of quantum field in the 
vacuum state 1 0(770)) is 



Po = ^ J d^kivkivo)- 



(62) 



This quantity is time dependent, but, considering the 
problem of particle oscillations in the present time, since 
the time scale of mixing phenomena are much smaller 
than the cosmological time scale, we can neglect the par- 
ticle creation in the vacuum due to the gravitational field 
and we may consider p ~ constant at the present time 
and then we renormalize in usual way. 

For further analysis of the vacuum structure in a 
curved background see refs.[l7l|. 



C. Mixed fields in expanding universe 

The boson mixing relations in FRW space-time are 
generalized as 

Xa(x,77) = xi(x,77) cos6'-|-X2(x,77) sin 61 
Xs(x,77) = -xi(x,77) sin6'-|-X2(x,77) cose* (63) 
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We now proceed in a similar way to what has been 
done in Ref. 10] for bosons in flat space-time and recast 
Eqs.|(nSl i^to the form: 



Xa{x,i]) = Gg ^(77) Geif]) 
Xs(x,?7) = Gg^{7]) X2(x,ry) Gg{i]) 



(64) 



and similar ones for 7ryi(x, 77), 7rs(x, 77). Gg{rj) denotes 
the operator which implements the mixing transforma- 
tions (jnSJ: 



Ge{v)^e^p[0Sirj)], 



(65) 



with 



S{v) = ~* y c^^x (^7ri(x,?7)x2(x,?7) 

- 7r2(x,77)xi(x,?7)y (66) 
We have, explicitly 

S{v) = J d'k(u^{r^)al,a^,2-V,r{v)a.u,ia^,2 

+ Vkiv) alj^ .al ^ ~ Ukiv) ak,ial^^) (67) 



where Ukir]) and Vk{i]) are Bogoliubov coefficients given 
by 

Ukiv) ^ ^ [vkAvMv) - v'kAvHM] ' (68) 

Vkiv) = -iWk,i{v)vk,2{v) ~ v'f,^2iv)vk,i{v)] ■ (69) 

and satisfy the relation 

\Uk\'-\Vk\' = 1. (70) 

Similar results can be obtained in the case of fermion 
mixing. 



For this space-time, the general solution of Ea. (|53() is 
Vkiv) = VWil \AJnim) + BY„{k\7j\)\ , (74) 



with n 



^ — j^, A and B constants, and J„(/c|77|) and 
y„(fc|77|) Bessel functions. 

In the de Sitter space-time, a suitable vacuum state 
is the Bunch-Davies vacuum, defined as the Minkowski 
vacuum in the early time limit (77 — )■ —00): 



Vkiv) 



1 



77 —> —00. 



(75) 



In this case, the mode functions are 



Vkiv) 



7r|77| 
2 



Jn{k\v\) - iYn{k\ri\) 



(76) 



with 77 = y'l - ^. 

Assuming for the present time, i.e. the actual time 
relative to the observer. 



t = 



77 



1 

H 



(77) 



the mixing transformations in the fiat space-time 
are good approximation of those in FRW space-time 
Eas. H63|l . since the time scale of mixing phenomena are 
much smaller than the cosmological time scale. The re- 
sults obtained in Section III describe the contribution to 
the value of the cosmological constant given by the neu- 
trino mixing at present time. 



V. CONCLUSIONS 



D. The de Sitter space-time 

The de Sitter space-time is characterized by a scale 
factor of the form 



a{t) 



(71) 



where H — a/a > Q is the Hubble constant. The con- 
formal time 77, the scale factor 0(77) and the effective fre- 
quency are given by 



77 



-Ht 



H 



0(77) 



1 

1 



77- 



(72) 
(73) 



The conformal time rj ranges from — cx) to when the 
proper time t goes from —00 to 00. The origin of 77 is 
chosen so that the infinite future corresponds to 7y = 0. 



The neutrino mixing is a possible source for the cos- 
mological constant. Indeed, the non-perturbative vac- 
uum structure associated with neutrino mixing leads to 
a non-zero contribution to the value of the cosmological 
constant The value of A consistent with its accepted 
upper bound is found by using the natural scale of the 
neutrino mixing cut-off. The origin of the contribution 
here discussed is completely different from that of the or- 
dinary contribution to the vacuum zero energy of a mas- 
sive spinor field: it is the mixing phenomenon which pro- 
vides the vacuum energy contribution discussed in this 
paper. 
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